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ABSTRACT 

Solar filaments exhibit a range of eruptive-like dynamic activity, ranging from the full or partial eruption of the 
filament mass and surrounding magnetic structure as a coronal mass ejection to a fully confined or failed eruption. 
On 201 1 June 7, a dramatic partial eruption of a filament was observed by multiple instruments on board the Solar 
Dynamics Observatory ( SDO ) and Solar-Terrestrial Relations Observatory. One of the interesting aspects of this 
event is the response of the solar atmosphere as non-escaping material falls inward under the influence of gravity. 
The impact sites show clear evidence of brightening in the observed extreme ultraviolet wavelengths due to energy 
release. Two plausible physical mechanisms for explaining the brightening are considered: heating of the plasma 
due to the kinetic energy of impacting material compressing the plasma, or reconnection between the magnetic field 
of low-lying loops and the field carried by the impacting material. By analyzing the emission of the brightenings 
in several SDO/ Atmospheric Imaging Assembly wavelengths, and comparing the kinetic energy of the impacting 
material (7.6 x 10 26 -5.8 x 10 27 erg) to the radiative energy (~1.9 x 10 25 -2.5 x 10 26 erg), we find the dominant 
mechanism of energy release involved in the observed brightening is plasma compression. 
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1. INTRODUCTION 

Solar filaments (called prominences when observed on the 
solar limb) exhibit a range of eruptive behavior, including 
dramatic activation with the filament mass remaining confined 
to the low corona (e.g., Ji et al. 2003; Alexander et al. 2006), 
the eruption of part of the observed filament structure (Zhou 
et al. 2006), and the almost complete eruption of all of the 
filament mass (Plunkett et al. 2000). The most common type of 
eruption is partial eruption, where prominence mass is observed 
falling back to the solar surface (see Gilbert et al. 2007 for a 
thorough discussion on different types of eruptions in the context 
of kinking motions). Prominence eruptions (full and partial) are 
often associated with coronal mass ejections and may play an 
important role in their initiation. 

Partial eruptions are particularly interesting because of what 
occurs after the eruptive part of the material escapes; the 
remaining supporting magnetic structure and prominence mass 
may relax or return to a lower altitude after the magnetic field 
reconfigures. This leads to the following question: what happens 
to the returning mass? Observationally, partial eruptions are 
most obvious on the limb, where pieces of material can be seen 
falling along apparent magnetic field lines in the plane of the 
sky (Gilbert et al. 2001), but are occasionally observed on the 
solar disk (Pevtsov 2002; Liu et al. 2012). Another phenomenon 
associated with filament eruptions that has historically been 
linked to returning material is the two-ribbon flare. One early 
explanation for two-ribbon flares is provided by Hyder (1967), 
in which the author concludes that two-ribbon flares are due 
to a chromospheric flare-like brightening mechanism. More 
recently, Shibasaki (2002) proposes that plasma falling from 
high altitudes after filament eruptions can convert potential 
energy into thermal energy, providing the source of energy and 
mass supply for long-duration solar flare events in the long- 
lasting decay phase. 

In the present work, we are interested in the physics behind 
the interaction of falling material and the solar atmosphere upon 


impact, which have the observational signature of brightening. 
Other phenomena with a similar observational signature are se- 
quential chromospheric brightenings (SCBs; Balasubramaniam 
et al. 2005; Pevtsov et al. 2007; Kirk et al. 2012). Kirk et al. 
(2012) propose that SCBs are caused by the chromospheric im- 
pact of accelerated plasma along newly reconfigured magnetic 
field lines. Further research will reveal whether the SCBs are 
related to the process described in this Letter, or if they are a 
fully independent phenomenon. 

The break-up of returning plasma after the 2011 June 7 
filament eruption (see Figure 1) and the fluid instabilities 
associated with the falling material were described recently by 
Innes et al. (2012). Williams et al. (2013) estimated the mass 
of the rapidly falling prominence material using high cadence 
extreme ultraviolet (EUV) images from the Solar Dynamics 
Observatory (SDO)/ Atmospheric Imaging Assembly (AIA), 
while Reale et al. (2013) used a combination of observations 
and simulations to investigate solar surface brightening due to 
the falling material. This event was also shown by Li et al. 
(2012) to trigger a globally propagating EUV wave. 

The focus of the present Letter is the interaction between 
falling material from the dramatic partial filament eruption and 
the solar atmosphere as the material returns to the surface. We 
address the physical mechanisms that are potentially responsible 
for the observed brightenings, with the objective of determining 
which mechanism is dominant. In this context, we utilize imag- 
ing observations from both SDO/ AIA and Solar -Terrestrial 
Relations Observatory (Sr£R£0)/Extreme-Ultraviolet Imager 
(EUVI). In doing so, we address the following question: are 
the observed EUV brightenings caused by the prominence ma- 
terial dissipating its kinetic energy in the chromosphere via 
collisions (compression), with the compressively heated plasma 
dissipated primarily by heat conduction, or are they the result of 
reconnection occurring between magnetic field lines involved in 
the impact (reconnection)? The type of reconnection postulated 
here is between the magnetic field carried by (or frozen into) 
the falling plasma and the ambient magnetic structure. Both 
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Figure 1 . SDO/PdK 193 A and STEREO -A /EUVI 195 A images with overlays showing the trajectories of four pieces of falling filament material following the 2011 
June 7 flare. The context images are shown at four time intervals, close to the respective impact times of each piece. In each panel, the track corresponding to impact 
region 1 is denoted by diamonds, region 2 by triangles, region 3 by squares, and region 4 by crosses. 


processes (i.e., compression and reconnection) are plausible but 
have different emission signatures. 

2. DATA 

Prominences are commonly observed above the solar limb in 
H (X = 6563 A) or Hei (10830 A) and Hen (10830 A) emis- 
sion. When seen projected against the solar disk, prominences 
appear as dark features in chromospheric lines — such as H 
and He i — due to absorption and also appear above the solar 
limb and against the disk in EUV coronal lines, such as Fe xn 
(195 A). For 195 A, the observed radiation arises from hydro- 
gen and helium continuum absorption. The filament eruption on 
2011 June 7, which occurred in NOAA AR 11226 (S22W55) 
was observed in multiple wavelengths by AIA (Femen et al. 
2012) onboard SDO (Schwer et al. 2002; Pesnell et al. 2012), 
which takes full-disk images in 10 (E)UV channels at Of! 6 spa- 
tial resolution and high temporal cadence of 12 s. The fila- 
ment material returning to the surface appeared in absorption in 
the EUV lines. It was also observed by EUVI (Wuelser et al. 
2004; Howard et al. 2008) imaging package onboard STEREO 
(Kaiser et al. 2008). EUVI provides observations in four pass- 
bands, namely, 171 A (Feix), 195 A (Fexn), 284 A (Fexv), 
and 304 A (Hen). The 195 A (Fexn) data was used for this 
study. 

3. ANALYSIS AND RESULTS 

The two mechanisms, compression and reconnection, are 
both viable explanations for the observed brightening observed 
in the SDO / AIA passbands upon impact of in-falling filament 
material. Clearly, the details of the magnetic topology of the 
falling material are not well-known; therefore, we use the term 
reconnection generically, indicating plasma heating by magnetic 
energy release. To determine which mechanism is responsible 
or dominant, we calculate the kinetic energy of the falling 
material and compare it with the energy associated with the 
observed emission. If the energies are comparable, evidence 
that the material is dissipating its kinetic energy via collisions, 
then the more likely mechanism is compression. If the energy 
release in the emission is larger than the kinetic energy, then 


reconnection is most likely playing an important role (i.e., the 
energy release is too large to be explained by compression). 
Allowing for uncertainties in the energy estimate, we impose a 
requirement that the energy release must be at least an order of 
magnitude larger than the kinetic energy to favor reconnection. 
This is reasonable and accounts for large uncertainties in the 
measurement methods. 

Although there is visual evidence that a large amount of 
prominence material is in-falling and interacting with the 
solar atmosphere, we restrict our analysis to five of the most 
accurate measurements and estimates of energy, observed from 
four distinct pieces of impacting prominence material. The 
trajectories of the impacting material are illustrated in Figure 1 
(see Figure 3 for precise impact areas). 

3.1. Kinetic Energy 

To determine kinetic energy, the velocity and mass of the 
prominence material impacting the chromosphere must be 
found. 

3.1.1. Velocities 

Line-of- sight velocity of the falling prominence material may 
differ substantially from the true velocity. Also, the presence 
of the magnetic field in the background corona may affect 
the trajectory of the material, which can deviate from purely 
ballistic trajectory, and the resulting kinetic energy estimate — 
proportional to v 2 — may differ by an order of magnitude de- 
pending on the line of sight. To obtain an accurate measure of 
the true velocities of the falling material (Figure 2), triangula- 
tion measurements were performed using observations from the 
STEREO-A spacecraft, combined with observations from SDO, 
utilizing the IDL routine in SolarSoft, ssc_measure (Thomp- 
son 2009). This enabled reconstruction of the three-dimensional 
path of the falling prominence material. 

As Figure 2 illustrates, the velocities of material prior to 
impacting the corona range from approximately 150 km s -1 
to over 300 km s _1 . We use these velocities in conjunction 
with mass measurements of the pieces of material to obtain an 
estimate of the kinetic energy associated with the impact. 
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Figure 2. Height vs. time plots showing the paths of tracked prominence 
material associated with each impact site (1-4; see Figures 1 and 3), ob- 
tained via triangulation measurements from combining SDO / AIA 193 A and 
STEREO-A/ EUVI 195 A images. The radial velocity for each piece of material 
was calculated from the final two height measurements. 
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3.1.2. Mass 

Gilbert et al. (2005) developed a technique for deriving 
prominence mass by observing how much coronal radiation 
in the Fexn (195 A) spectral line is absorbed by prominence 
material. In the present work, we apply this method, which 
allows us to consider the effects of both foreground and 
background radiation in our calculations, to the pieces of 
prominence falling from the 2011 June 7 eruption to obtain a 
measure of prominence density. This method also accounts for 
“volume blocking” or the amount of coronal radiation that would 
be present where the prominence is located. Volume blocking is 
described in detail by Heinzel et al. (2008) and Labrosse et al. 
( 2010 ). 

The calculation of prominence column mass density 
(g cm -2 ) along the line of sight requires the determination of 


the extinction factor (a measure of how much coronal radiation 
is being absorbed as it travels through a prominence). If o is the 
mean absorption cross-section for radiation passing through a 
prominence, the extinction factor for radiation traveling in the 
direction s over a distance l is 


a = e ti nads . 


( 1 ) 


where n is the total prominence number density. If o is uniform 
throughout the prominence, and if we define the column density 
by 



( 2 ) 


then from Equations (1) and (2) we have 


N = -a~ l In a. (3) 

To obtain a, we measure the intensity in the SDO / AIA 193 A 
images in the region of falling prominence material and the 
region just outside each side of that same prominence material 
(to interpolate background intensity behind the prominence 
material). Due to the nature of the equations, the above- 
mentioned intensity measurements are taken in two adjacent 
regions characterized by very different background intensities 
(denoting these regions by the superscripts L and D for light 
and dark), which occurs at time 07:15 UT as a large piece of 
material is seen crossing the solar limb (see Gilbert et al. 2005, 
2006 for a complete derivation). 

Assuming that the foreground radiation in the light region is 
simply related to that in the dark region by a proportionality 
factor that we can specify: 

//=£//. ( 4 ) 

With some manipulation of the equations found in the full 
derivation, we obtain expressions for the extinction factor and 
the foreground radiation: 


E - PE 

O' = —R T- 

A? - PE 


( 5 ) 
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Figure 3. Top row: AIA 193 A images of the chosen brightenings associated with prominence material impact. The white contour in each panel shows the defined 
impact area. Bottom row: AIA flux data from each impact region, as defined in the top panel. All six optically thin EUV wavelengths are shown. 
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Table 1 

Summary of Impact Region Properties 


Impact 

Start Time 
(UT) 

End Time 
(UT) 

Solar Location 
(arcsec.) 

Mean Impact Area 
(10 17 cm 2 ) 

Radial Velocity of 
Material (km s -1 ) 

K. E. 
(10 27 erg) 

f^rad 

(10 25 erg) 

la 

lb 

2 

3 

4 

08:07:50 

08:08:50 

07:26:30 

07:03:00 

07:51:30 

08:09:00 

08:11:12 

07:31:00 

07:07:00 

07:56:00 

(407, -75) 
(393, -65) 
(445, -385) 
(825, 25) 
(475, 100) 

0.9 

3.2 

6.6 

12.8 

2.8 

-200 ± 15 
-200 ± 15 
-183 ±36 
-146 ± 13 
-305 ± 37 

0 76 +0 - 72 
U - /O -0.45 

2.70+ 2 j 7 5 6 9 

4.66^ 2 

c 7 c+6.15 

c 4Q+6.6I 

D -^-3.46 

1 0+1-3 
i,y -0.7 

8-42x0 

24.72 7 4 7 1 

4 4+1-4 
4>4 -0.8 
0 o+0.9 
Z,8 -0.5 


Applying this technique to the SDO/ AIA 193 A data at 
07:15:45 UT and assigning uncertainties to the various com- 
ponents, we find 

« = 0 . 028 ^ 0*0206 ’ ( 6 ) 

and 

taa = -(3.561+V2w). (7) 

Representing the fractional hydrogen and helium abundances 
(by number) by /h and /n e (where /h + fn& ~ 1), and the 
H and He ionization fractions by Xu = ^h/(^h + ^h + ) and 
*He = ^He + /(^He + ^He + X we Can write 

= /h(1 — ^h)<Th + /ke(l — •^He)° r He + /ke*HeO r He+ (8) 

where a H , (XHe> and a He + are the photoionization cross-sections 
(for 193 A radiation) for H, He, and He + . The mean prominence 
mass corresponding to the mean cross-section in Equation (8) 
is m = /ke^He + /h^H = (4 /He + /h)^H- 
Taking cr H = 7.69 x 10 -2 ° cm 2 , cr He = 1.54 x 10 -18 cm 2 , 
and crHe+ = 1.01 x 10 -18 cm 2 (T. A. Kucera & H. R. Gilbert 
2013, in preparation; Ready & Kilcrease 2000), we find from 
Equation (8) that a = (1.83 ±0.4) x 10 -19 cm 2 (see Section 5 in 
Gilbert et al. 2005 for the range of ionization states considered), 
and from Equations (3) and (7) the column density is 

N = (l.94+'|f 5 |) x 10 19 cm -2 . (9) 

Williams et al. (2013) use multi- wavelength SDO imaging 
data to apply polychromatic and monochromatic methods in 
estimating column hydrogen densities of pieces of falling 
filament material from the same event. Their estimates are 
larger (AhIO 20 cm -2 ) for a targeted piece of material, which 
they describe as exceptional, but are consistent with our results 
for the thinner threads of material (i.e., lower limit values in 
the range Ah ^ 10 18 -10 19 cm -2 ), which describe the material 
analyzed in the current work. 

If the mean mass per particle in the prominence is m, then 
the column mass density p is given by fi = Nm , and we can 
integrate fi over the prominence area (as seen by the observer) 
to obtain the prominence mass M p : 



In the present work, we evaluate a at one point in the 
prominence, and we must estimate an effective prominence area 
A e ff, such that 


7 / 


M p = (4/ He + fn)mn / / Nda « (4/ He + /h)^h x AA eff . 


(ID 

Based on the brightening areas of the impacting material 
(see Table 1 below), and assuming the blobs of material being 


measured have the same width and depth dimension, we estimate 
A q ff. We expect this to be an upper limit on the area due to 
the possible spreading of emission after initial impact. Since 
mu = 1.67 x 10 -24 g, it follows from Equations (9) and (11(7)) 
that the prominence mass range for the tracked pieces of material 
is 

M p = 1.82 x 10 12 -9.42 x 10 13 g. (12) 

Given the range of velocities of the impacting pieces of 
material in Section 3.1.1 and the estimate of the mass of those 
pieces in Equation (12), we obtain the kinetic energy K.E. 
of the impacting material. We find that K.E. is in the range 
7.58 x 10 26 -5.75 x 10 27 erg. The full K.E. estimates and 
uncertainties are listed in Table 1 . 

3.2. Observations of Plasma Impacts and Estimation 
of the Radiated Energy 

In order to estimate the radiated energy from the observed 
brightenings at the impact sites, we use the observed AIA fluxes 
from the brightenings to estimate the differential emission mea- 
sure (DEM) of the associated plasma volume. This procedure is 
based on the approach taken by Aschwanden et al. (2013). The 
properties of this DEM are used to estimate the total radiated 
energy per second dL^/dt. 

In this analysis, five of the most clearly observed impacts 
were selected for study. Figure 3 shows the location of each 
of these impacts, two of which (la and lb) are associated with 
the same piece of infalling prominence material. Figure 3 also 
shows the light curves of each of the labeled impacts in each 
optically thin AIA wavelength band. 

For each of these impacts, we use the recorded flux in each 
AIA channel in order to estimate the DEM distribution as a 
function of time for the plasma associated with the impact 
brightening. As described above, this is achieved via forward 
modeling, where we choose a distribution of the form, 

dEM f'logr-iogr^ 

~dT~ = EM 0 exp ( J, (13) 

i.e., a Gaussian emission measure distribution with peak temper- 
ature T c and width a , as utilized by, e.g., Aschwanden & Boemer 
(2011), Aschwanden & Schrijver (201 1), and Aschwanden et al. 
(2013). 

The temperature response functions of the AIA channels are 
the source of significant uncertainty and remain the subject 
of active study (e.g., Aschwanden et al. 2013), particularly 
the 94 A and 131 A channels at low temperatures. In order to 
account for this, we include a 25% uncertainty in the measured 
AIA flux at each wavelength due to instrument response, as 
suggested by Boerner et al. (2012) and Guennou et al. (2012). 
This is combined in quadrature with the statistical uncertainty 
associated with the AIA flux measurements. 
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Figure 4. Examples of the fitting of differential emission measure distributions to the AIA fluxes from each impact region. The top panels show the ratio of the fitted 
to the observed flux in each channel; the bottom panels show the best-fit DEM as a function of temperature. 


The best fit to the observed flux is achieved at each time 
interval via a search over the parameter space given by the 
variables EMq, T c , and a using the x 2 test. Figure 4 shows an 
example best-fit for each impact, found in each case near the 
time of peak EUV emission. 

A reasonably constrained DEM allows us to estimate the total 
radiated energy from the emitting plasma. The radiative loss rate 
may be written (e.g., Aschwanden 2005), 

— — = / EM(T) x A (T) dT erg s" 1 , (14) 

dt J Tl 


where A (T) represents the radiative loss function and EM{T) 
is the emission measure multiplied by the emitting area A, 
and hence is in units of cm -3 . For each emitting region, A 
is estimated from the AIA observations of the brightening 
extent (see Figure 3). The radiative loss function A (T) has been 
previously estimated by many authors, and varies primarily as 
a function of temperature (White et al. 2005). For this work, we 
obtain an estimate of A (T) from the CHIANTI database (Dere 
et al. 1997; Landi et al. 2012), choosing appropriate coronal 
abundances. 
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The last step toward estimating the total radiated energy is to 
integrate over the time duration of the event, hence, 

Uad = f ' dL 'f n dt. (15) 

Jt B dt 

For each impact, the dashed lines in Figure 3 denote the 
start and end times of the integration. Following this proce- 
dure, we find that the total radiated energy of these events is 
10 25 -10 26 erg. The full results are summarized in Table 1. 

To estimate the uncertainty associated with the total radiated 
energy estimates, a number of Monte-Carlo simulations were 
performed on the calculation of dL v ^/dt. For each simulation 
run, 1000 best-fit DEMs were generated based on simulated 
AIA input fluxes. These were sampled from a distribution 
generated using AIA flux measurements of the impacts and 
their uncertainties. The resulting uncertainties on the radiated 
energy estimates are shown in Table 1 . 

4. DISCUSSION 

The energy release in emission is one to two orders of 
magnitude smaller than the calculated kinetic energy. Given 
our requirement that the energy release must be at least an 
order of magnitude larger than the kinetic energy to favor 
reconnection, our results show a clear indication that the 
compression mechanism dominates in producing the observed 
brightenings. Moreover, we note an observational signature 
supporting compression over reconnection: all wavelengths in 
the SDO analysis respond at the same time and have similar 
decay times. This is not expected in the case of reconnection 
since the various lines would show different response times 
as the plasma cools after being heated. One might argue that 
secondary pulsations seen in the AIA data are a signature 
of flaring, but we note that in several cases the impacting 
prominence material is extended into multiple pieces that impact 
the same area, causing an expected temporal variation in the 
deposition of energy. 

Although these results strongly support that one mechanism 
is dominant over the other, both are likely occurring since the 
falling material undoubtedly carries frozen-in magnetic flux. Ex- 
amining more examples of falling prominence material would 
provide a more complete understanding of the relative impor- 
tance of the two mechanisms. Studying the temporal evolution 
of the emission and comparing the emitted energy to the impact- 
ing kinetic energy can provide insight on the efficiency of energy 
conversion and on the properties of the plasma, such as density, 
filling factors, thermal, and magnetic pressures. By comparing 
the emissivity and the evolution of these events to emission 
from flares in events not associated with falling material can 
help us understand the relative importance of the magnetic dis- 
sipation and the compression that occurs through the Lorentz 
force acting on the flaring plasma, ultimately improving our 
understanding of coronal heating. A three-dimensional MHD 
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model of an active region impacted by falling cool and dense 
prominence material needs to be developed in order to improve 
the understanding of these events. These important tasks are left 
for future studies. 
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